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Abstract — Many numerical simulations in quantum (bilin- 
ear) control use the monotonically convergent algorithms of 
Krotov (introduced by Tannor in [12]), Zhu & Rabitz ([11]) 
or the general form of Maday & Turinici ([13]). This paper 
presents an analysis of the limit set of controls provided by 
these algorithms and a proof of convergence in a particular 
case. 

I. INTRODUCTION 

The control of quantum phenomena is a tòpic that has been 
(and is still) a source of many interesting challenges not only 
to physics and chemistry but also to the mathematics and 
applied mathematics communities ([1], [2]). At the level of 
the experiments, làser control of complex molecular systems 
is becoming feasible, especially since the introduction ([3], 
[4]) of closed loop laboratory learning techniques and their 
successful implementation ([5], [6], [7], [8], [9], [10]). 
On the other hand, at the level of the numerical simulations, 
the introduction of the monotonically convergent algorithms 
by Zhu & Rabitz ([11]) that extend an algorithm due to 
Krotov ([12]) has allowed a considerable progress and made 
possible further investigations in this area. Recently, a gen- 
eral class of monotonically convergent algorithms has been 
proposed ([13]) and a relevant time discretization has been 
developed ([14]). 

However, no general analysis to explain in depth the conver- 
gence of these algorithms is available to date. In an attempt 
to fill this gap, this paper presents some results on the set of 
the controls provided by monotonic algorithms. 
Note that, among others, this question was raised in [15], 
but a wrong statement about the Cauchy character of the 
sequence is made that makes the proof not working as stated; 
the proof is more involved as explained in what follows. 
The balance of the paper is as follows: the necessary back- 
ground and definitions of the quantum control settings are 
given in Section II, properties of the monotonic sequences 
are presented in Section III, followed by the properties of the 
limit set in Section IV. Further results in a particular case are 
given in Section V and concluding remarks are presented in 
Section VI. 

II. QUANTUM OPTIMAL CONTROL AND 
MONOTONIC SCHEMES 

A. Cost functional and Euler-Lagrange Equations 

Consider a quantum system prepared in an initial state ipo 
and whose dynamics is characterized by its internal Hamil- 
tonian H. By assumption this Hamiltonian does not give 
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rise to an appropriate evolution and an external interaction 
is introduced in order to obtain the desired final property. 
This interaction is taken here as an elèctric field with time- 
dependent ampli tude e(í) that influences the system through 
a time-independent dipole moment operator /i. The new 
Hamiltonian H — fie(t) gives rise to the equations (we work 
in atòmic units i.e. h = 1): 



i—ip{x, t) 



dt 

ip(x, t 



0) 



Hip(x,t) 



fj,(x)e(t)ip(x, t) 



where we denote by x the relevant spatial coordinates. These 
equations hold on SI = R. N but for numerical tests we will 
consider that x belongs to an interval íl = [0, L] and that 
ip(0,t) — ip(L,t) = 0, for a large enough real number 
L and any t in R. This approach is justified by physical 
reasons since wave functions are generally localized in a 
space interval. 

The optimal control framework is then used to find a 
suitable evolution of e(t). The goal that the final state i/j(T) 
has prescribed properties is expressed by the introduction of 
a cost functional J to be maximized. This cost functional also 
includes a contribution that penalizes undesirable effects. 
One simple example of such a cost functional is: 



J{e) = {^{T)\0\^{T))-a 



2 {t)dt, 



(D 



where a > is a parameter (it may also depend on time cf. 
[16], [17]) and O is an observable operator that encodes the 
goal: the larger the value (tp(T)\0\ip(T)) is, the better the 
control objectives are met (here and in what follows we use 
the convention that for any functions / and g and any op- 
erator F: (f\F\g) — J f(x)Fg(x)dx. Note that, in general, 
achieving the maximal possible value of (ifj(T)\0\ip(T)) is at 
the price of a large làser influence J Q T e 2 (t)dt ; the optimum 
evolution will therefore strike a balance between using a low 
làser fluence while simultaneously maximizing the desired 
observable. 

At the maximum of the cost functional J(e), the Euler- 
Lagrange critical point equations are satisfied ; a Standard 
way to write these equations is to use a Lagrange multiplier 
x(x,t) called adjoint state. The following critical point 



equations are thus obtained ([11]): 

i^(x,t) = {H - e{t)fj)^{x,t) 
i/)(x,t = 0) = Vo(z) 

i&X&t) = (H - e{t)n) X (x,t) 



and assume that it has been proven that \\e 
M, ||er fe — 1 II 2 < M. Since we also have: 



2 < 



X (x,t = T) = (hl>(x,T) 
ae(t) = -Im(x(í)|M|V(í)). 



(2) 

(3) 
(4) 



From now on, ip(t) and \(t) w iU represent two functions of 
the Hilbert space L 2 (Í7; C) for almost all t in [0,T]. 

B. Definition of the monotonic schemes 

Efficient strategies for solving in practice the critical 
point equations (2)-(4) are represented by the monotonically 
convergent algorithms ([11], [12], [13]) that are guaranteed 
to improve the cost functional J at each iteration. In the 
formulation proposed in [13], the monotonic algorithms are 
described by the resolution of the following equations at step 
k: 



i^ k {t) = {H-e k (t) l· i)xP k {t) 
#(z,í = 0) = Vo(z) 



(5) 



e k {t) = (1 - 5)t~ l {t) - -Im( x fe - 1 (í)|Ml^(í)}(6) 



il tX k (t) = (H-ï k {t)n) X k {t) 
X k (x,t = T) = 04> k {x,T) 



(7) 



È*(t) = (1 - V )e"{t) - -±lm( x k (1)1^(1)). (8) 
a 

where S and rj are two real parameters. 

The most important property of this algorithm is given in 
the following theorem ([13]): 

Theorem 1: Suppose O is a self-adjoint positive semi- 
definite operator. Then, for any rj,ó € [0, 2] the algorithm 
given in Eqns. (5)-(6) converges monotonically in the sense 
that: 

J{e k+1 ) > J(e k ). (9) 

III. PROPERTIES OF THE SEQUENCE (e k ) k , {è k ) k 

We first prové that {e k )k and ie k )k defined in (5) and in 
(6) are bounded. We then prové that every weakly convergent 
subsequence is strongly convergent. In the following, ||.| 
represents the norm of L 2 (Í7, C), whereas |j.|| 2 represents 
the norm of L 2 ( [0, T] ; R). The scalar product in L 2 ( [0, T] ; M) 
will be denoted by < ., . >. 

A. Bound for the sequences 

We suppose from now on that O and \x are bounded 
operators and we denote by ||0||», IHI* tne i r norms. 

Theorem 2: There exists M > such that, for all k > 0, 
the solutions e k , s of (5-8) verify: 

Ví £ [0,T], |e fc+1 (í)| < M, |e fe+1 (í)| < M. 

Proof: Define M by: 
M=max(|| £ || 2 , ||g°|| 2) max(l, * ^BMl), 



\\s% < ii - m+ -ii 1 1 ► iimix'-Hm^m h, 

the Cauchy-Schwartz inequality yields: 

Kx'-Hm^m < Hx^WII·IIM^W)!! 

< IH|..||x fc - 1 (í)l|.||^ fc (í)ll· 

We then use the following equalities and bounds on state and 
adjoint state: 

Ví,||^(í)|| = 1, 
WJIx^WII - \\x k - 1 {T)\\ = \\OW> k - 1 (T))\\ 
< ||0||..||V fc - 1 (T)|| = ||0||. 

to obtain the estimate: 

\\e%<\l-S\M + S m *- M * . 

a 

If 6 < 1, then the definition (10) yields ll°ll*JIHI. < M and 
then: ||e fc || 2 < |1 - S\M + 5M = M , and if 5 > 1 then 
_s_\\ojï^Mk < M and in this case: ||£ fe || 2 < |1 - S\M + 

={S- 1)M + (2 - S)M = M. 
A similar proof leads to the same estimate for e*. ■ 

B. Weak convergence of subsequences 

1 ) Extraction of a weakly convergent subsequence: Be- 
cause e k is bounded in the Hilbert space L 2 ([0, T]; M), there 
exists a weakly convergent subsequence that will be denoted 
by (e kn ) n . Let e be the weak limit associated to (e fe ")„. 

2) Límits of (e k " +1 - e fe ")„ and (è*" +1 - £*»)„: The 
sequence J(e k ) is bounded since |J(e fe )| < ||0|| + M. It 
has also been proven ([13]) that: 

J{e k+1 )-J(s k )=(ip k+1 (T) - ip k (T)\0\ip k+1 (T) - ip k (T)) 
+ j\l-l)(e k +\t)-è k {t)fdt 

+ [ T (--l)(è k (t)-e k (t)fdt, 
Jo V 

which gives after summation: 

N-l 

J(e N )-J(e°)=J2(<P k+1 (T) ^ k (T)\0\i> k+1 (T) ^(T)) 
o 

N-l 



í à-l^e^^-t^fdt 
Jo S ^ 

Ç T 9 N-l 

/ (--l)£(g*(í)- £ fc (t)) 2 dí. 
JO V n 



2-5' 



(10) 



Thus the series J^o^ \\s k+1 and ||è* -£ k \\í 

converge and we deduce that: 

lim||e fe " +1 -e fe "|| 2 = lim ||g*" +1 - l k ™ || 2 = 0. (11) 

n n 

Similar results hold when 5 = 0,rj ^ and 5 = 1, 77 7^ 0. 
Remark: Such properties do not guarantee the convergence 



of the sequences. For example, the sequence (u n ) n defined 
by u n = sin(log(n+l)) verifies E^K+i-* 1 ») 2 < 
however it does not converge. 

3) Weak convergence of (e kn+p ) n and (è kn+p ) n : Let è 
be a test function in L 2 ([0, T]; R). From: 

< i, e kn+1 >=< è, e kn+1 - e kn > + < s, e kn >, 

one can easily prové that (e kn+1 ) n weakly converges to e 
too. By the same way, (e kn+p ) n and (s kn+p ) n also weakly 
converges to e. 

C. Strong convergence of (sk n ) n 

1) Strong convergence of (ip kn ) n , (x fe " _1 )n and (x )n-' 
Since we have proven that {ek n )n and (sk n )n weakly 
converge in L 2 ([0, T]; R), hence in L 1 ([0, T]; R), we can 
use theorem 3.6 of ([18]), which implies that ip kn strongly 
converges in C([0, T]; L 2 (íí, C)) to the state ip associated to 
e. One can also easily adapt the proof of this theorem in order 
to obtain that {x kn ~ 1 )n an d (x fe ")« also strongly converge 
in C([0, T]; L 2 (fí, C)) to the adjoint state associated to s and 
4>{T). 

2) Strong convergence of (e kn ) n : The strong con- 
vergences of (ip kn ) n , {x kn ~ 1 )n an d (x kn )n implies 
the strong convergence of (— lm(x kn ~ 1 \l· L \' ] P kn ))n an d 
(fIm(x fc "|MlV' fe "})r l in C([0,T];R). According to the def- 
initions (5) and (6), we can now write (s kn ) n as follows: 

e k " +1 = (l-6)(l- V ) e k "+u n , 

S v ' 

A 

where (u n ) n strongly converges. Let e denote a positive real 
number. Since |A| < 1, there exists an integer jo sucn that 
A- 70 1 < e. Let us write then: 

io-l k„-j 
j=0 j=0 

The second term of (12) is a finite sum of strongly convergent 
contributions and its third term can be bounded above by 
e ll u ll2ï3x> which ends the proof of the strong convergence 
of (e fe ")„. The strong limit is necessarily e. Passing to the 
limit in (5), we conclude that e is a critical point of J. 
A similar proof can be done to prové that (s kn ) n strongly 
converges to e. 

It thus has been proven that every weakly convergent subse- 
quence of (e k )k strongly converge in C([0, T]; R) to a critical 
point of J. 

IV. PROPERTIES OF THE LIMIT SET 

This section is devoted to the study of the limit set of 
the sequence (e k )k which will be denoted by A. From now 
on we will suppose that this set contains at least two points 1 . 

'Of course A depends on the initial guess s° for the definition of the 
sequence (e k )k (i.e. A = A(e°) even though we shall skip this dependency 
in all what follows). 



A. First properties 

From theorem 2, one deduces that A C B(0, M), where 
M is defined in (10) and 5(0, M) stands for the ball of 
radius M of L 2 ([0, T];R). According to the definition of A, 
the results of the latest section prové that A is a subset of the 
set of critical points of J. Finally, thanks to the monotonic 
property (9) of the sequence (e k )k, we deduce that J is 
constant on A. 

B. Compactness 

Let us now prové a first topological property. 

Lemma 1: The set A is compact. 

Proof: Let (ff^Jn denote a sequence of A. We can 
associate to this sequence a subsequence (e fc ")„ of (e k )k 
such that || — e fc "||2 < ^. According to the previous 
results we can extract from (e fe ")„ a strongly convergent 
subsequence (e fc ™')„/. Let e* denote the limit of this latter. 
Thus, the sequence (e^Jn' strongly converges to e* that is 
a point of A. ■ 

C. e-Strings in A 

Consider a general mètric space (E,d), (x,y) E E 2 and 
e a positive real number. We call e-string between x and y 
a finite sequence z\,...,Zn of point of E such that: 

• Z\ = x, 

• z N = y, 

. Vfc G [Í,N- 1], d{z k +\z k ) < e. 
Then the set A has the following topological property. 

Lemma 2: For any (e^,^) e A 2 and any e > 0, there 
èxits an e-string in A between and e'^. 

Proof: As a compact set, there exist N open balls 
of radius f covering A. By the definition of A and (11), 
there exists an infinity of K > for which Ik — 
£oo, £ K , e K+N ( K \ e'^ is an e-string. From Ik, one can 
then build another e-string l' K = e K ' 1 ,e K ' 2 , ...,s K ' N °. In- 
deed, if N a > N(K), define l' K by: 

/' _ 7^ F K+N(K) K+N (K) 

s v ' 

No-N(K) terms 

and if No > N(K), one can remove N(K) — Nq terms of 
Ik while keeping the e-string properties. 
For each i, 1 < i < N , let us extract from {e k ' 1 )k 
a strongly convergent subsequence of {e k )k- The limits 
obtained are an e-string in A. ■ 

D. Connexity 

The previous result leads to the following theorem. 

Theorem 3: The set A is connex. 

Proof: Suppose there exist two closed subsets of A, 
denoted by A\ and A 2 , such that A = A\ U A 2 and Ai D 
A 2 = 0. Because of the existence of e-strings for every e, we 
deduce that the distance between A\ and A 2 is equal to 0. 
Since A is compact, this is in contradiction with A\ n^b = 0- 



E. Summary 

It has been proven that the limit points of a sequence 
obtained by a monotonic scheme are a compact and connex 
set of critical points of J. Note that if this set is reduced 
to one point, the compactness of the sequence implies its 
convergence. 

V. VARIATIONAL ANALYSIS AND PARTICULAR 
CASE 

Let us focus now on the scheme obtained for 5 = 1 and 
i] = 0, which corresponds to the Krotov formulation (as 
in [12]). We will estimate the variations of ip and \ with 
respect to e. The results obtained will enable us to prové the 
convergence for large vàlues of the parameter a. 
The above defined set A is still considered to contain at least 
two points. 

A. Estimates 

Let e and e' be two points of A, ip and the correspond- 
ing states given by (2) and x an d x' me corresponding ad- 
joint states solution of (3). Consider (2) written in integrated 
form, for ip and tp': 

Jo 

V>'(í) - e- lHt ^ + f e- lH{t - s h'{s)i l· L^{s)ds. 
Jo 

Let us introduce the notations Sip(t) = ip(t)—tp'(t), Sx(t) = 
x(t) — x'{t) an d Se(t) = e{t) — e'(t), we then have: 

6ip(t) = í e- iH ^6e(s)inip(s)ds 
Jo 

+ [ e- iH ( t -'h , (a)iii6il)(8)d8. (13) 
Jo 

Since the operator e~ lHt is unitary, we deduce that: 

\\íe- iH ^Se{s)i^(s)ds\\<\\^\UT\\d£\\ 2 , 
Jo 

1 1 f e- iH{ - t -^e , {s)i l· i5^{s)ds\ \<M\ |, / 1 \8^{a) \ \ds, 
Jo Jo 

where M has been defined in (10). From Gronwall's lemma 
applied to (13), we obtain: 

\\Sm\\ < |N|.Te T ll"ll· M ||fe||i ) (14) 

where ||.||i represents the norm of L 1 ([0, T]; M). A similar 
computation for the adjoint state leads to: 

||<5x(í)ll<HO||*||MlUT(l + e T IWI* M )e T IWI· M ||fc|| 1 . 

(15) 



B. Convergence 

Since e and e' are critical points of J, the two following 
equalities hold: 

ae(t) = -Im( X (í)|M|V'(í)}, 
as'(t) = -Mx'WImI/W)- 
The difference of these two equalities yields: 

aSe(t) = -Tm{(6 X \n\m) + (xWm))- 
From (14,15) we have: 

a||&||i<||0|UIHI^ 2 (l + e T l^l· M )e 2T ^H* M ||fc|| 1 . 

(16) 

Thus we get the following result: 

Theorem 4: The monotonic scheme defined by (5)-(8), 
6 = 1, r] = strongly converges in L 2 ([0,T];R) under the 
assumption that: 

a>\\0\UMlT 2 (í + e T ^' M )e 2T ^ M . 

Proof: Suppose that the monotonic scheme does not 
converge, then there exists at least two distinct points e and 
e'. Using the above notations, the equation (16) holds in this 
case. Since Se ^ 0, we reach a contradiction. ■ 

VI. CONCLUSION 

It has been proven that the sequences provided by mono- 
tonic schemes are compact and that the set of their limit 
points is compact and connex. It has been shown that this set 
reduces to one point (i.e. the algorithm strongly converges) 
for a large làser fluence penalty parameter a. We refer the 
reader to [19] for a more detailed presentation of this tòpic. 
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